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We demonstrate a correlation between the transmissivity of sandstones to visible light and their fluid 
permeability. By means of experiments carried out on sandstone and sand packs, we show that the light 
transmission increases when the contrast between the indices of refraction of the matrix material and 
pore-filling fluid is reduced, and that, for a given fluid, it is greatest for rocks with the largest grains. We 
present qualitative and analytical arguments, and results of a numerical simulation, which account for 
these correlations. The simulations also yield the space and time distribution of the transmitted light. 
Finally, we show that light transmission provides a method of mapping the permeability of sandstone 
without fluid flow, and at high spatial resolutions. Extension of these findings to rocks other than sand- 


stones is discussed. 


I. INTRODUCTION 


The nondestructive study of rock interiors employs 
many types of energy—acoustic, nuclear, x-ray, and 
microwave—but generally not visible light. Light is not 
used because it is common experience that visible light 
cannot penetrate a rock. On the other hand, the indivi- 
dual grains of many consolidated rocks are composed of 
nonabsorbing material. Thus, if a rock composed of such 
grains is filled with a pore fluid having a sufficiently simi- 
lar index of refraction, one expects that light will indeed 
penetrate such a rock. 

In this paper, we present the results of an experimental 
study of sandstones using visible light. We find not only 
that light can sometimes be transmitted through sand- 
stone, but also that the transmission coefficient is corre- 
lated with an important aspect of the pore structure of 
the rock—its permeability to fluid flow. In particular, 
we find that large permeability is correlated with low at- 
tenuation of visible light. In view of this correlation, we 
suggest a new method for mapping the permeability of a 
porous rock without fluid flow and at a spatial resolution 
not previously realized. 

The permeability k of porous rock is defined by the 
Darcy relation: v= —(k/1)Vp, where v is the fluid ve- 
locity, 7 is the fluid viscosity, and p is the hydrostatic 
pressure. A large permeability is of central importance 
to the extraction of oil.!'* In a rock formation, k typically 
has considerable spatial variation. In studies of such for- 
mations, one is interested in k at many length scales, 
from millimeters to kilometers. Permeability is usually 
measured by causing a fluid (typically oil, water, or air) to 
flow through a portion of rock under the influence of a 
pressure gradient, the rock sample being contained in a 
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suitable fashion.>~°> These measurements are laborious 
and slow, involve boundary conditions that are difficult to 
interpret,° and usually achieve resolutions no finer than a 
few cm. Attempts have been made to characterize rock 
permeability on a smaller scale’~° and to measure direc- 
tional permeability in anisotropic rock.'°!! 

Composite media have been extensively studied by op- 
tical probes in recent years, especially in the short- 
wavelength limit.!2~!? Optical probes of permeability 
may not at first appear to be a suitable alternative to con- 
ventional methods, even if light does penetrate a rock. 
For instance, fluid flow is confined to the pore space, 
while, if the rock grains and pore fluid are both transpar- 
ent, light penetrates both rock and pore space. However, 
rock permeability is strongly affected by typical pore di- 
mensions: fine pores are typically correlated with low per- 
meabilities due to large viscous drag forces created by 
high surface-to-volume ratios. Similarly, for a given pore 
fluid, and with grain dimensions d large compared to the 
wavelength A, light will scatter more strongly inside a 
fine-grained (i.e., fine-pored) rock than in a large-grained 
one because of the greater number of scattering surfaces. 
Light is therefore expected to be more strongly attenuat- 
ed in fine-grained composities. This argument suggests a 
simple correlation between permeability and light propa- 
gation: rocks with low fiuid permeability should also 
tend to attenuate light strongly. 

In the next section, we describe light propagation ex- 
periments carried out on two types of rocklike compos- 
ites. In Secs. III and IV, we present a theoretical analysis 
of these experiments, using both qualitative arguments 
and numerical calculations to account for the experimen- 
tal trends. Next, we describe optical experiments on nat- 
urally occurring rocks (sandstones) and show that there is 
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indeed a strong correlation between optical transmission 
and permeability. We conclude with a demonstration of 
high-resolution optical mapping of the permeability of a 
sandstone sample. 


I, EXPERIMENTS 
ON UNCONSOLIDATED COMPOSITES 


To deal with relatively controlled geometries, we have 
initially studied two types of unconsolidated composites, 
chosen to resemble natural sandstones: dense packings of 
spherical glass beads, and packs of rounded, but irregu- 
larly shaped, quartz sand grains. Microscopic examina- 
tion revealed that the glass beads were clear, with smooth 
spherical surfaces, while the quartz grains were transpar- 
ent with a slight brownish tinge and multifaceted sur- 
faces”? that were optically “rough” (i.e. had surface irre- 
gularities of a scale much smaller than d).”’ The refrac- 
tive index 7 of glass is 1.52 at visible wavelengths, while 
that of quartz is 1.55.2? Three size-sorted glass bead 
packs were used. These had size ranges 100-110 ym (re- 
ported as d=105 pm), 170-180 wm (175 pm), and 
250—300 um (275 wm). Four sand packs were used, with 
size ranges 53-74 um (64 wm), 74-105 um (90 pm), 
88—125 um (107 wm), and 177-250 wm (214 wm). The 
word “size” for the glass beads refers to the diameter of 
the spherical beads. For the sand grains, the sizes quoted 
here are defined by the mesh sizes of the sieves used to 
generate the various batches of sand. 
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FIG. 1. Experimental setup: (a) sample compartment and (b) 
measurement geometry. 
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The glass bead or sand pack was contained in a com- 
partment with glass walls. A 10-mw He-Ne laser was 
used as a light source (red light, wavelength 632.8 nm, 
randomly polarized; beam diameter 0.68 mm). The 
transmitted optical power was measured by a photodetec- 
tor calibrated for the 632.8-nm wavelength. The area of 
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FIG. 2. Optical attenuation of a laser beam of wavelength 
632.8 nm through a slab of thickness 2.6 mm of glass bead pack 
of spherical beads (a) and of sand grain packs comprised of 
rounded quartz grains (b). Results are shown for several sizes of 


glass beads or quartz grains saturated with three fluids: (i) air, 


(11) water, and (iii) toluene. Continuous curves are the results of 


numerical simulations. 
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the photocell, placed flush against the glass wall of the 
sample compartment, was approximately 1X1 cm? (Fig. 
1). To eliminate the effects of reflection at various inter- 
faces, two sample thickness L (2.6 and 5.2 mm) were 
studied. The ratio of the two levels P, and P, of 
transmitted power, expressed in decibels, is reported as 
the attenuation A, L.e., 


A=10.0log,(P,/P,) dB . (1) 


The attenuation includes both the loss due to beam 
scattering, as well as any actual absorption in the constit- 
uents. In our samples, however, the measured loss is al- 
most entirely due to the former cause: true absorption is 
very small. 

Measurements were carried out on packs saturated 
with three fluids: (i) air, (ii) brine with a resistivity of 0.2 
Q, m, and (iii) a fluid consisting of 93% toluene and 7% 
paraffins of refractive index similar to toluene. In case 
(ii), a Separate experiment showed that the results of the 
optical experiment were the same with either brine or dis- 
tilled water. The various saturating fluids have indices of 
refraction 1 ;, = 1.00, 1 ater = 1-33, and Mioiuene = 1-49." 

The results of the attenuation measurements are 
presented in Fig. 2. The experimental uncertainty of the 
measurements is estimated to be less than 1 dB. These 
results suggest three main conclusions: (i) the attenuation 
is dependent on the size parameter d of the beads or 
grains, smaller sizes being correlated with larger attenua- 
tion, (ii) the attenuation, for a given d, decreases as the 
refractive indices of the matrix and pore fluid become 
more similar, and (iii) for a given d, and given pore fluid, 
the attenuation is generally much larger for the quartz 
grains than for the glass beads. 

In order to have an independent estimate of the porosi- 
ties of the sand and bead packs, the dc resistivity of each 
brine-saturated bead and sand pack was measured under 
identical conditions. These resistivities were found to be 
nearly equal for all the glass or sand packs considered, 
suggesting that the porosities @ (volume fraction occu- 
pied by pore space) were also approximately equal. Thus, 
not surprisingly, the results of Fig. 2 cannot be explained 
as a dependence of transmission on porosity. 


Ill. QUALITATIVE THEORETICAL ARGUMENTS 


We now present a simple mean-free-path argument 
that qualitatively accounts for the general features 
(i)—(iii) listed above. The transmission coefficient T 
through a slab (defined as the transmitted power divided 
by the incident power) depends on, among other parame- 
ters, the slab thickness L, the indices of refraction n, and 
n, of the rock matrix and the pore fluid, the typical linear 
dimensions d of the pore or rock grain, and the wave- 
length A of the incident electromagnetic wave in free 
space. 

To develop a qualitative estimate of the attenuation, we 
assume that we are in the short-wavelength limit A <d, 
which is appropriate for the packs and the frequency con- 
sidered here. We also assume that the pores have sur- 
faces that are flat on a scale of A. A ray traveling 


through this medium will, on average, encounter a 
pore/grain interface approximately every distance d, be- 
ing refracted according to Snell’s law: n,sin0,;=n,sin@,, 
where 6, and @, are the angles made by the ray in medi- 
um 1 and in medium 2 with the surface normal. The 
justification for considering only the refracted ray 1s 
given in Sec. IV. Writing 6,=06,1+A0, and expanding 
Snell’s law to first order in AO, we obtain 


Ad=(n, —~Hy tan, /n., ’ 


which is valid to first order in An/n,,, where 
An=n,—n), and n,, is the volume-averaged index of re- 
fraction. The quality n,, here could equally well be taken 
as n, or n, to the same order in An /n,,. 

We now use this formula to obtain an estimate of the 
number of scatterings required before a given ray has lost 
the memory of its original direction—that is, has been 
scattered through an angle ~7/2 rad. Since @, is a ran- 
dom number depending on the (random) orientation of 
the surface to the incoming ray, tan@, is also a random 
number, of order unity. Thus, to first order in An /n,,, 
each scattering produces a change of angle of order 
A@~An/n,,. Since the scattering angle is equally likely 
to be positive or negative, the root-mean-square change 
of angle (A@), produced after N scatterings will be given 
by a random-walk-like formula, i.e., (A@)y=AnVN / 
ny. The rms number of scattering events N, required to 
scatter a given ray through an angle 7/2 will be of order 
N.=[n,,7(2An)~']*. Hence, the “elastic photon mean 
free path” /.,, defined as the distance a photon (i.e., a 
light ray) will travel before it is scattered through an an- 
gle 7/2, is approximately 
2 


(2) 











The experimental results (1) and (ii) summarized in the 
preceding section can be qualitatively understood from 
this relation. It states that /,, (a) increases with pore or 
grain size, and (b) increases with decreasing |An|. Since 
the transmission 7 through the slab increases with in- 
creasing /,,, it follows that, for a fixed slab thickness, T 
should increase with increasing pore size and with de- 
creasing |An|. 

Our understanding of the experimental result (i11) is 
somewhat more speculative. The glass bead pack consists 
of beads of uniform size and shape, with no surface 
roughness such as characterizes the sand grains. In ap- 
plying relation (2) for /., to grains with rough surfaces, 
the grain diameter d should more appropriately be re- 
placed by the correlation length € characterizing the sur- 
face roughness of the grain—that is, the typical distance 
a photon must travel, on average, before encountering a 
pore-grain interface. It is shown in the Appendix that a 
better estimate for /,, is 

2 
la ee ee ‘ (3) 
((S€)°) In(2k7é*) 
where k=27V/€,,/A is the wave number, and de is the 
difference between the dielectric constants of the two 
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constituents. Thus, except for logarithmic corrections, 
/,, increases linearly with the correlation length €. For 
the glass bead packs, & is probably of the order of the 
bead diameter, or larger, since the surfaces of the beads 
are smooth. On the other hand, for the quartz grains, & 
is probably smaller than the grain diameter because of 
variations in grain shape and surface roughness. Hence, 
this simple formula suggests that, for glass beads and 
quartz grains of the same nominal diameters, the 
transmission through the quartz packs should be smaller, 
as is observed. 


IV. NUMERICAL SIMULATION 


We next consider a simulation of propagation in the 
geometrical optics limit (A <<d), with zero absorption, 
through a slab of thickness L and infinite extent in the 
two transverse directions. A ray enters the sample at 
normal incidence, and propagates a random distance x, at 
which point it is assumed to encounter a pore-grain inter- 
face (represented as a randomly oriented planar surface). 
We follow either the refracted or the reflected ray with 
probability R or 1—R, where R is the power reflection 
coefficient, computed by averaging over polarizations 
perpendicular and parallel to the plane of incidence: 


R=4([r,|7+|r,/7) ’ (4) 


r, and r, being the usual Fresnel amplitudes for reflection 
of a ray polarized parallel and perpendicular to the sur- 
face.*? Since usually R <<1, we generally follow the re- 
fracted ray. 

After scattering from one surface, the randomly chosen 
(reflected or refracted) ray travels another random dis- 
tance and meets another randomly oriented surface. The 
reflected or refracted ray from this surface is again select- 
ed at random as described above, and the process contin- 
ues until the ray exits from the front or the back of the 
sample. If it reaches the detector, it is counted and its 
energy, location, and time of arrival are noted. Each ray 
begins with an energy © normalized to unity, which is re- 
duced upon each scattering by a factor of R for a 
reflected ray or 1—R for a refracted ray. The transmis- 
sion coefficient T is calculated as 


T=fO6gy » (5) 


where f is the fraction of followed rays that reach the 
detector, and &,, is the average energy of a ray that does 
reach the detector.”4 

Figure 2(b) shows our results for a sample of quartz 
grains saturated with three different fluids corresponding 
to those in the experiment. In each case, the average dis- 
tance X,,,i;, traveled by a ray within a grain between 
scatterings is taken as 1.5d, where d is the typical grain 
size defined in Sec. II. We choose a value slightly larger 
than d, since a ray can travel from one grain to another 
without passing through the intervening pore space, i.e., 
by passing through a grain-grain contact area. The dis- 
tance between scatterings is taken as a random variable 
uniformly distributed between 0 and 3d. We choose the 
typical distance x,,. traveled in pore space to be 


0.69X1.5d, uniformly distributed between O and 
1.38X1.5d. We take x joe <X grain because the estimated 
porosity ¢ of the experimental samples is less than 0.5, so 
that rays typically travel smaller distances in the pore 
space. The factor of 0.69 is obtained from the assump- 
tion that 6=0.25, which gives 


X pore /X grain =16/(1—$)]'*° = 0.69 : 


We typically consider 10° rays, which is sufficient to en- 
sure the convergence of T to within a few percent. We 
deduce the attenuation, as in the measurements, by calcu- 
lating the transmissions T, and T, through two slabs of 
thickness 2.6 and 5.2 mm. The calculated value of 
10log,9( 7, /T, ) is quoted as the attenuation (in dB).”° 

The results of the simulation show reasonable agree- 
ment with the data for quartz sands over the entire range 
of grain sizes and for the three types of pore fluids. In 
particular, the simulation correctly reproduces the in- 
crease in attenuation with decreasing grain size, and with 
increasing contrast between the indices of refraction of 
pore fluid and grains. 

Figure 3 shows the calculated spatial distributions of 
the transmitted light for quartz saturated with water for 
several particle sizes. In each case, we have plotted I(r), 
defined as the transmitted intensity at z=L=2.6 mm 
and at a radial distance r from the axis of the laser beam, 
given an incident beam normally incident on the slab at 
r=Q, z=0 of a cylindrical coordinate system. The figure 
is normalized so that f I(r)(27r )dr is proportional to the 
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FIG. 3. Calculated transmitted intensity I(r) for 2.6-mm- 
thick quartz sandpacks of various average grain diameters sa- 
turated with water, plotted as a function of the radial distance r 
from the center of the detector (beam axis). The various intensi- 
ties are normalized so that for a given particle size, the areal in- 
tegral f I(r )2tr dr equals the transmission coefficient T. 
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total transmission coefficient JT for each grain diameter. 
The half width of the distribution is broadest for the 
smallest particles, where the rays have been most strong- 
ly scattered before emerging from the slab. We have ob- 
tained similar graphs for the spatial intensity distribu- 
tions for quartz saturated with other fluids: in general, 
for a given slab thickness and particle size, we find, as ex- 
pected, the spatial distribution is widest when there is the 
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greatest contrast between the indices of refraction of 
quartz and saturating fluid. 

In Fig. 4 we plot the distribution of travel times of the 
rays for a slab of thickness 2.6 mm containing quartz 
grains of various sizes and saturated with three different 
fluids. The figure plots P(7), the transmitted power (in- 
tegrated over the plane z=ZL) at time 7 after the laser 
beam is incident at the front of the slab, normalized so 
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FIG. 4. Calculated transmitted power P(r) plotted as a function of arrival time 7, integrated over the surface at z=L, for quartz 
sandpacks of various average particle sizes saturated with toluene (a), water (b), and air (c). Normalization is such that for each sa- 
turating fluid, the area under a curve equals the transmission coefficient T. The absolute normalization is arbitrary but is the same 
for each of the three sets of curves. 
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that f P(r)dr=T, the transmission coefficient. If the 


transmission were ballistic, P(+) would be a 6 function at 
T=L/(c/n,,), where c/n,, is the speed of light in the 
medium. The width of the distribution once again is seen 
to increase with increasing scattering within the slab, i.e., 
with decreasing grain size or increasing index-of- 
refraction contrast.”° 

The temporal and spatial distributions of arriving 
pulses suggest other ways of analyzing the transmission 
through the slab, which may, like the transmission itself, 
yield information about the permeability, grain size, pore 
fluid, and other geometrical properties of the medium. 
Thus an experimental arrangement that could measure 
these distributions might be of interest in analyzing such 
properties of the composite. 


V. APPLICATION 
TO NATURALLY OCCURRING ROCKS 


As shown in the previous sections, for a given saturat- 
ing fluid, the transmission is typically largest for quartz 
sand packs with the largest grains—the same materials 
for which the fluid permeability is expected to be largest. 
Figure 5(a) shows that this expected correlation is also 
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FIG. 5. (a) Transmitted optical power plotted against per- 
meability for five samples of Boise sandstone. Experimental ar- 
rangement is shown in the inset. The sample disks, 5 mm thick, 
are saturated with water and kept immersed in water up to a 
level 2 mm above the upper surfaces of the samples. The depar- 
ture from a monotonic correlation exhibited by point B is be- 
lieved to be due to deviations in thickness from its nominal 
value. (b) Transmitted optical power plotted against porosity, 
for the same samples as in Fig. 5(a). No correlation is observed. 


observed in Boise sandstone, a relatively clean (clay-free) 
sandstone in which the spatial variation of fluid permea- 
bility arises from a spatial variation in grain size. The air 
permeabilities were measured on inch-long, inch- 
diameter cylindrical plugs, from which 0.2-in.-thick disks 
were cut. This thickness is small enough to permit a 
measurable optical transmission through the sample, 
while at the same time being sufficiently large to be 
characterized by the permeability of the parent plug. 
Figure 5(a) displays the measured optical power transmit- 
ted through water-saturated samples, plotted against the 
measured permeability of the samples with respect to the 
flow of air. There is a clear correlation—large per- 
meabilities correspond to large optical transmissivities. 
Figure 5(b) shows the optical transmission plotted as a 
function of the sample porosity, showing that these quan- 
tities are uncorrelated. 

Given a theory connecting fluid permeability to light 
transmission, or a calibration diagram such as Fig. S(a) 
for a given type of rock, we can hope to substitute, for 
some purposes, the optical measurement for the direct 
measurement of the permeability. The optical measure- 
ment can be made with high spatial resolution. Since the 
laser beam spreads out into an approximately conical 
volume as it propagates through the sample, this resolu- 
tion is thickness dependent, the best resolution being ob- 
tained for thin samples. 

To illustrate the potential use of this method, we have 
examined a slab of Berea sandstone 10.2 cm square and 
initially 2.54 cm thick. The sample is homogeneous on 
large scale but has a stratified structure on a centimeter 
to millimeter scale with the layers perpendicular to the 
slab. These layers differ in grain size. A conventional 
permeability map for the slab was first obtained using a 
standard air minipermeameter.? The measurements were 
made on grids 1.27 cm square and have a resolution of 
about 2 cm. 

The slab was then sliced down to a thickness of 5.1 mm 
SO as to permit appreciable laser transmission upon satu- 
ration with water, while still remaining thick enough to 
retain bulk permeability characteristics as measured by 
the minipermeameter. The same face of the slab was 
presented to the laser beam as to the minipermeameter 
nozzle. The optical image is shown in Fig. 6(a), taken on 
grids 1.27 mm square. The experimental arrangement is 
similar to that shown in the inset of Fig. 5(a), with the 
sample translated by micrometer stages to obtain areal 
coverage. The layering is clearly evident. Using a cali- 
bration diagram such as Fig. 5(a), one could, in principle, 
convert the optical transmission to permeability values. 
This suggests a means of utilizing the optical method to 
develop an automated permeability mapping system. 

In Fig. 6(b) we show the data of Fig. 6(a) replotted 
after calculating the harmonic mean over every 100 grid 
squares so that the new grids correspond to those of the 
miniperimeter measurements. (Arithmetic and geometric 
means produce similar diagrams.) Data for the perimeter 
grid points are omitted to avoid anomalous values due to 
edge effects. Figure 6(c) is the permeability map obtained 
using the minipermeameter. Within the uncertainties of 
the measurements, the resemblance between Figs. 6(b) 
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FIG. 6. (a) Optical transmission map of a slab of Berea sandstone. (b) Data of (a) replotted upon averaging over every 100 grid 
squares. (c) Permeability map for the same slab obtained using an air minipermeameter. 


and 6(c) is satisfactory; in particular, the presence of the 
layers is evident in both sets of measurements. This ex- 
ample suggests how optical and conventional permeabili- 
ty measurements can be combined to test schemes for 
analyzing the permeability of porous rocks at various 
length scales. 

The same method could be applied, in principle, to 
other types of rock, such as the carbonates, comprised of 
limestones and dolomites. In such rocks, unlike clean 
sandstones, the grains are opaque. Any transmitted light 
is therefore largely confined to the pore space. Thus we 
still expect a grain-size-dependent transmission, but prob- 
ably not a strong dependence of transmission on the re- 
fractive index of the pore fluid. Our preliminary experi- 
mental studies show that a detectable optical transmis- 
sion still occurs through limestone, though the attenua- 
tion is much greater than in sandstone. The carbonate 
rocks are therefore best studied using more sensitive 
detectors or higher-power lasers. Such rocks may exhibit 
an even stronger correlation between fluid permeability 
and light transmission than the sandstones, since in the 
carbonates, both fluid and light flow are confined to the 
pore space. 


VI. CONCLUDING REMARKS 


We have shown that a correlation exists between opti- 
cal transmission through artificial composites and natural 
rocks and the grain size or permeability of such materi- 
als. We have also presented a simple theoretical analysis 
of the results. Many other potential applications of the 
methods are conceivable, both for natural rocks and for 


media other than natural rocks. Both the experiments 
and the theoretical interpretations presented here are rel- 
atively simple: further refinements of both are necessary 
to provide a detailed quantitative understanding. For ex- 
ample, the interpretation of experiments is complicated 
by pore-filling and pore-lining substances such as clay. 
Another complication, not discussed here, is the effect of 
topology.*’ In measurements such as described here, in 
which the contrast between the indices of refraction of 
grain and pore fluid is not too great, light is transmitted 
through both media and the influence of pore topology 
on light transmission is probably not critica]. In cases of 
greater contrast, this topology may be more important, 
possibly leading to an even stronger correlation between 
optical transmission and. fluid permeability. Finally, we 
note that light of different wavelengths may elicit other 
types of information, especially if the relevant indices of 
refraction are frequency dependent. 
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APPENDIX: CALCULATION OF /,, 
FOR A RANDOM COMPOSITE 
IN THE SHORT-WAVELENGTH LIMIT 


We want to compute an elastic mean free path for light 
scattering from a random composite in the short- 
wavelength limit. Thus, we imagine a volume V of ran- 
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dom composite, embedded in a homogeneous medium of _eikr 
dielectric constant €) (which may be a function of fre- D,(r)= —~ fe, Eo) » (A3) 


quency). Writing the dielectric function at point r as = 
e(r)=e,(r)+Se(r), we have a linear relation between where k is the scattered wave vector, k =oV/ €,)/c (elastic 


D(r) and E(r) (suppressing the frequency dependence): scattering), and € and €) are unit polarization vectors 
parallel to the electric field vectors of the scattered and 
D(r)=€9E(r)+ de(r)E(r) . (Al) — incident waves, respectively. In the Born approximation, 


Sn the scattering amplitude f(€,é,)) takes the form 
We now assume that the composite is subject to an in- ae (€,€p) 


cident plane wave with electric field £(2,8,)=— tix . f de(r’) Dire d3r'. (Ad) 
E,(r,t)=E,expl(iky:r—iot) , (A2) ae “0 

where ke Ve je: Where werea homonencous medi: The corresponding differential cross section is 

um of dielectric constant €, filling all space, the electric de. S| £(E;eq)I-) 

displacement would be simply D)p>=€,)Ep. Because of the dD = Ip? : (A5) 

presence of the composite, the total displacement is : 

D=D,+D,, where D, is an additional scattering term. where (---) denotes an ensemble average over the 


In the Born approximation, and at distances far from the configurations of the random function de(r’). Substitut- 
composite in comparison to its linear dimensions, it can ing Eq. (A4) into Eq. (A5), and averaging over polariza- 


be shown from Maxwell’s equations that tions of the incident radiation, one obtains 
be Ee Se ee es ee 
4 
do — oul 24 1 + cos7@) f ( Se(r’ )Se(r’’) Ve Tig es ytd apn (A6) 
dQ 32m’ 


ee 
where q=k—ko, 0 is the angle between k and ky, and the __ that the mean free path /., is given by 














integral runs over the volume V of the composite. 1 1 Ae 
To calculate the differential scattering cross section, 7. = ag |) osha 0 (A9) 
one must know the form of the correlation function el 
(Se(r' )de(r’’)). If one assumes an exponential correla- _ k4*((8e)*) rr mip 
tion function of the form a 22 fl, 
(8e(r)5e(r’)) = (Se)? exp —— (A7) where 
3+ 8k7E?+ 8k 4E4 
where & is an appropriate correlation length (assumed to L= 16k 8¢8 In(1+4k 76") 
be much smaller than the dimensions of the composite), 
one obtains _— (34+14k7E + 16k 4E") ALD 
: 4k®E(1+4K7E?) 
da _ Vk*((8e)*) E*(1+cos76) 
dQ Are? [1+2k2E%(1—cos6)]? | (A8) In the long-wavelength limit (k*é*<<1), [=4 and 


1/1, is proportional to k*, as expected for Rayleigh 
scattering. In the short-wavelength limit (k?é?>>1), we 


The elastic mean free path is the distance traveled by a find 


photon before it loses memory of its original direction. 
By analogy with the corresponding calculation in the 1 ((Se)*) In( 2k *é?) 

theory of electronic mean free paths,”® this is obtained by cae E , (A12) 
weighting the differential-scattering cross section by a . 

factor (1—cos@), which favors large angles. We then find with is equivalent to Eq. (3). 


1A. I. Levorsen, Geology of Petroleum (Freeman, San Francisco, (Academic, New York, in press). 
1967). 7B. F. Swanson, J. Pet. Technol. 33, 2498 (1981). 

2J. Bear, Dynamics of Fluids in Porous Media (American E]- 8A. J. Katz and A. J. Thompson, Phys. Rev. B 34, 8179 (1986). 
sevier, New York, 1972). °M. A. Chandler, D. J. Goggin, and L. W. Lake, J. Sediment 

3R. E. Collins, Fluid Flow Through Porous Materials (Reinhold, Petrol. 59, 613 (1989). 
New York, 1961). 104. E. Scheidegger, Geofis. Pura Appl. 28, 75 (1956). 

4A. E. Scheidegger, The Physics of Flow Through Porous Media 11F.M. Auzerais, D. V. Ellis, S. M. Luthi, E. B. Dussan, and B. 
(University of Toronto Press, Toronto, 1974). J. Pinoteau (unpublished). 

°F, A. L. Dullien, Porous Media: Fluid Transport and Pore 12K. M. Yoo, Feng Liu, and R. R. Alfano, Phys. Rev. Lett. 64, 
Structure (Academic, New York, 1979). 2647 (1990). 


6Y. Bernabé, in Fault Mechanics and Transport Properties of 13... F. Gate, J. Phys. D 4, 1049 (1971). 
Rocks: The Brace Volume, edited by T. F. Wong and B. Evans 14A. Ishimaru, Y. Kuga, R. L. T. Cheung, and K. Shimizu, J. 


204 B. R. DE, I. H. H. ZABEL, D. STROUD, AND M. A. NELSON 45 


Opt. Soc. Am. 73, 131 (1983). 

154. Ishimaru and Y. Kuga, J. Opt. Soc. Am. 72, 1317 (1982). 

164. Z. Genack, Phys. Rev. Lett. 58, 2043 (1987). 

17D, J. Pine, D. A. Weitz, P. M. Chaikin, and E. Herbolzhei- 
mer, Phys. Rev. Lett. 60, 1134 (1988). 

18M. P. Van Albada, J. F. de Boer, and A. Lagendijk, Phys. 
Rev. Lett. 64, 2787 (1990). 

19S. Feng, L. Golubovic, and Y. C. Zhang, Phys. Rev. Lett. 65, 
1028 (1990). 

270W. A. Deer, R. A. Howie, and J. Zussman, Rock-forming 
Minerals (Wiley, New York, 1963), Vol. 4. 

214. Ishimaru, Wave Propagation and Scattering in Random 
Media (Academic, New York, 1978), Vol. 2. 

22 Handbook of Chemistry and Physics, 53rd ed., edited by R. C. 
Weast (Chemical Rubber Co., Cleveland, 1972), pp. E-208 
and E-209. Crystalline quartz is birefringent, with principal 
indices of refraction 1.55 and 1.54 in the visible. We have ar- 
bitrarily assumed isotropy with n = 1.55 in our modeling. 

233, D. Jackson, Classical Electrodynamics (Wiley, New York, 
1975), Chap. 7. 

24R. Jullien and R. Botet, J. Phys. (Paris) 50, 1983 (1989) have 
used a procedure similar to ours to compute the scattering 
cross section of fractal aggregates, in the limit when the parti- 
cles making up the aggregates are large compared to the 
wavelength. The Monte Carlo criterion used by these authors 
for following the reflected or refracted wave is slightly 


different from ours but should lead to the same weighting of 
transmitted and reflected waves in the limit of a sufficiently 
large number of incident rays. 

25Thus, we do not assume that the transmitted intensity falls off 
exponentially with thickness. In fact, even for a perfectly 
nonabsorbing medium in the diffusion limit, the transmission 
does not fall off exponentially. See, e.g., P. W. Anderson, Phi- 
los. Mag. B 52, 505 (1985); D. J. Durian, D. A. Weitz, and D. 
J. Pine, Science 252, 686 (1991); or Ishimaru, Wave Propaga- 
tion and Scattering in Random Media (Ref. 21). 

26For a sufficiently thick sample, the width of pulses such as 
shown in Fig. 4 is determined by the diffusion coefficient for 
light, D=v_gl,,/3, where vg is a suitable velocity of energy 
transport [see M. P. Van Albada, B. A. Van Tiggelen, A. 
Lagendijk, and A. Tip, Phys. Rev. Lett. 66, 3132 (1991)]. 
Such pulse profiles have been calculated by a number of au- 
thors by solving the diffusion equation. The tails of these dis- 
tributions can be analyzed to give the absorption coefficient of 
the medium. See, for example, G. H. Watson, P. A. Fleury, 
and S. L. McCall, Phys. Rev. Lett. 58, 945 (1987). 

27A. H. Thompson, A. J. Katz, and C. E. Krohn, Adv. Phys. 
36, 625 (1987). 

28See, for example, J. M. Ziman, Principles of the Theory of 
Solids (Cambridge University Press, Cambridge, 1964), pp. 
187f. 


